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Abstract: 

Higgs-boson production in association with W or Z bosons, pp — > WH/ZH + X, is 
the most promising discovery channel for a light Standard Model Higgs particle at the 
Fermilab Tevatron. We present the calculation of the electroweak 0(a) corrections to these 
processes. The corrections decrease the theoretical prediction by up to 5-10%, depending 
in detail on the Higgs-boson mass and the input-parameter scheme. We update the cross- 
section prediction for associated WH and ZH production at the Tevatron and at the 
LHC, including the next-to-leading order electroweak and QCD corrections, and study the 
theoretical uncertainties induced by factorization and renormalization scale dependences 
and by the parton distribution functions. 
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1 Introduction 



The search for Higgs particles pQ is one of the most important endeavours for future high- 
energy collider experiments. Direct searches at LEP have set a lower limit on the Standard 
Model (SM) Higgs-boson mass of M H > 114.4 GeV at the 95% confidence level (C.L.) 0. 
SM analyses of electroweak precision data, on the other hand, result in an upper limit of 
M H < 211 GeV at 95% C.L. The search for the Higgs boson continues at the upgraded 
proton-antiproton collider Tevatron [3] with a centre-of-mass (CM) energy of 1.96 TeV, 
followed in the near future by the proton-proton collider LHC [5J with 14 TeV CM energy. 
Various channels can be exploited at hadron colliders to search for a Higgs boson. At the 
Tevatron, Higgs-boson production in association with W or Z bosons, 

pp^WH + X and pp^ZH + X, (1.1) 

is the most promising discovery channel for a SM Higgs particle with a mass below about 
135 GeV, where decays into bb final states are dominant |4,. 

At leading order, the production of a Higgs boson in association with a vector boson, 
pp — > VH + X, (V — W, Z) proceeds through qq annihilation [U], 

qq' -> V* -> V + H. (1.2) 

The next-to-leading order (NLO) QCD corrections coincide with those for the Drell-Yan 
process and increase the cross section by about 30% jjj. Beyond NLO, the QCD corrections 
for VH production differ from those for the Drell-Yan process by contributions where the 
Higgs boson couples to a heavy fermion loop. The impact of these additional terms is, 
however, expected to be small in general p], and NNLO QCD corrections should not 
increase the VH cross section at the Tevatron significantly, similar to the Drell-Yan cross 
section j^j. As described in more detail in Section IH the renormalization and factorization 
scale dependence is reduced to about 10% at 0(a s ), while the uncertainty due to the 
parton luminosity is less than about 5%. At this level of accuracy, the electroweak 0(a) 
corrections become significant and need to be included to further improve the theoretical 
prediction. Moreover, the QCD uncertainties may be reduced by forming the ratios of the 
associated Higgs-production cross section with the corresponding Drell-Yan-like W- and Z- 
boson production channels, i.e. by inspecting a p p^vH+x I 'cpp-^v+x- In these ratios, higher- 
order electroweak effects should be significant. For the Drell-Yan-like W- and Z-boson 
production the electroweak corrections have been calculated in Refs. [IHIIII] and [T2] . 
respectively. 

In this paper we present the calculation of the electroweak 0(a) corrections to the 
processes pp/pp — > W + H + X and pp/pp — > ZH + X. 1 We update the cross-section pre- 
diction for associated WH and ZH production at the Tevatron and at the LHC, including 
the NLO electroweak and QCD corrections, and we quantify the residual theoretical un- 
certainty due to scale variation and the parton distribution functions. 

The paper is organized as follows. In Sect. |2] we outline the computation of the 0(a) 
electroweak corrections. The calculation of the hadronic cross section and the treatment of 

1 The electroweak 0(a) corrections to associated ZH production at e + e~ colliders have been presented 
in Ref. [T3]. 
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Figure 1: Lowest-order diagram for qq' — > V* — * VH (V = W, Z). 



the initial-state mass singularities are described in Sect. El In Sect. 0]we present numerical 
results for associated WH and ZH production at the Tevatron and at the LHC Our 
conclusions are given in Sect. El 



2 The parton cross section 

2.1 Conventions and lowest-order cross section 

We consider the parton process 

9(p ? ,T ? ) + g'(p g -<,r r ) -> V(p v ,X v )+H(p H ) [+7(M)]> (2.1) 

where V = W + , Z. The light up- and down-type quarks are denoted by q and q', where 
q = u, c and q' — d, s for W + H production and g = g' = u, d, s, c, b for ZH production. 
The variables within parentheses refer to the momenta and helicities of the respective 
particles. The Mandelstam variables are defined by 

S=(Pq+Pq') 2 , t=(Pq-Pv) 2 , U = {p ql - p V f , S V H = {.P H + Pvf ■ (2.2) 

Obviously, we have s = syu for the non-radiative process qq' — > VH. We neglect the 
fermion masses m q , rrigi whenever possible, i.e. we keep these masses only as regulators in 
the logarithmic mass singularities originating from collinear photon emission or exchange. 
As a consequence, the fermion helicities r q and Tg> are conserved in lowest order and in the 
virtual one-loop corrections, i.e. the matrix elements vanish unless r q = —r q i = r = ±1/2. 
For brevity the value of r is sometimes indicated by its sign. 

In lowest order only the Feynman diagram shown in Fig. Q contributes to the scattering 
amplitude, and the corresponding Born matrix element is given by 

2 T 

Ml = e f'^r e(p<r) tv&v) "M Pq ), (2.3) 

where ey(Xy) is the polarization vector of the boson V, v(p q /) and u{p q ) are the Dirac 
spinors of the quarks, and u± = |(1 ± 75) denote the chirality projectors. The coupling 
factors are given by 

r V ud r „ My/ 

9udw = ~7k — 0t ~> 9wwh = , 

9 qq z Q q + — — d T -, gzzH - — — , (2.4) 

c w c w s w c w s w 
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where Q q and I? = ±1/2 are the relative charge and the third component of the weak 
isospin of quark q, respectively. The weak mixing angle is fixed by the mass ratio Mw/M%, 
according to the on-shell condition sin 2 6 W = sfy = 1 — c 2 ^ = 1 — M^/M|. Note that 
the CKM matrix element for the ud transition, V u d, appears only as global factor \V u d\ 2 in 
the cross section for WH production, since corrections to flavour mixing are negligible in 
the considered process. This means that the CKM matrix is set to unity in the relative 
corrections and, in particular, that the parameter V u d need not be renormalized. The same 
procedure was already adopted for Drell-Yan-like W production [TUl lllj. 

The differential lowest-order cross section is easily obtained by squaring the lowest- 
order matrix element Ai^ of ()2.3|) . 

fda \ 1 1 A 1 / 2 (M 2 ,Ml,s) _ . WT|2 



V dfi / 12 64tt 2 § 
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a 2 



48 Mis 1 



9vvh ((4'vO 2 + (W) 2 ) AV2 ( M v> s) 



(t - Mp(u - Ml) + M 2 v s 

(s - Ml) 2 ' l2 - 5j 



where the explicit factor 1/12 results from the average over the quark spins and colours, 
and Cl is the solid angle of the vector boson V in the parton CM frame. The total parton 
cross section is given by 



(X 7T 

a (qq -> VH) = i^pp 9vvh {{^'vf + 



72M 2 s 2 

xAi, '' tf ^^> >w 'S)' 2M ^ (2 - 6) 

where A is the two-body phase space function A(a;, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 
The electromagnetic coupling a = e 2 /(4:7i) can be set to different values according to 
different input-parameter schemes. It can be directly identified with the fine-structure 
constant a(0) or the running electromagnetic coupling a(k 2 ) at a high-energy scale k. For 
instance, it is possible to make use of the value of a(M|) that is obtained by analyzing ^3] 
the experimental ratio R = cr(e + e~ — > hadrons)/(e + e~ — > n + n~). These choices are called 
a(0)-scheme and a (M|) -scheme, respectively, in the following. Another value for a can 
be deduced from the Fermi constant G^, yielding q;g m = V^G^Ml^sl^/n; this choice is 
referred to as G ^-scheme. The differences between these schemes will become apparent in 
the discussion of the corresponding 0(a) corrections. 



2.2 Virtual corrections 

2.2.1 One-loop diagrams and calculational framework 

The virtual corrections can be classified into self-energy, vertex, and box corrections. 
The generic contributions of the different vertex functions are shown in Figs. |2] and El 
Explicit results for the transverse parts of the WW, ZZ, and ^Z self-energies (in the 



3 




Figure 6: Diagrams for the corrections to the qqH vertices 



't Hooft-Feynman gauge) can, e.g., be found in Ref. [T3]. The diagrams for the gauge- 
boson-fermion vertex corrections are shown in Figs. HJ El and El The diagrams for the 
corrections to the WWH, ZZH, and jZH vertices can, e.g., be found in Figs. 8 and 9 
of Ref. [12] • The box diagrams are depicted in Figs. and |H1 where (p is the would-be 
Goldstone partner of the W boson. 

The actual calculation of the one-loop diagrams has been carried out in the 't Hooft- 
Feynman gauge using standard techniques. The Feynman graphs have been generated 
with FeynArts [T71[TB] and are evaluated in two completely independent ways, leading to 
two independent computer codes. The results of the two codes are in good numerical 
agreement (i.e. within approximately 12 digits for non-exceptional phase-space points). In 
both calculations ultraviolet divergences are regulated dimensionally and IR divergences 
with an infinitesimal photon mass m 7 and small quark masses. The renormalization is 
carried out in the on-shell renormalization scheme, as e.g. described in Ref. |15j . 

In the first calculation, the Feynman graphs are generated with FeynArts version 1.0 
|17j . With the help of Mathematica routines the amplitudes are expressed in terms of 
standard matrix elements, which contain the Dirac spinors and polarization vectors, and 
coefficients of tensor integrals. The tensor coefficients are numerically reduced to scalar 
integrals using the Passarino-Veltman algorithm |19j . The scalar integrals are evaluated 
using the methods and results of Refs. j201[TH] . 

The second calculation has been made using FeynArts version 3 jTHj for the diagram 
generation and FeynCalc version 4.1.0.3b [21] for the algebraic manipulations of the am- 
plitudes, including the Passarino-Veltman reduction to scalar integrals. The latter have 
been numerically evaluated using the Looptools package [22] version 2. 

2.2.2 Renormalization and input-parameter schemes 

Denoting the one-loop matrix element in 0(a) the squared matrix element reads 



\M T + M\\ 2 = \M T \ 2 + 2 Re{(M r YMl} + ... = (1 + 2 Re{^ irt }) |7W-| 2 + . . . . (2.7) 



Substituting the r.h.s. of this equation for |-Mq| 2 in (j2.5|) includes the virtual corrections 
to the differential parton cross section. The full one-loop corrections are too lengthy and 
untransparent to be reported completely. Instead we list the relevant counterterms, all 
of which lead to contributions to M.\ that are proportional to the lowest-order matrix 
element M- T ct = S^Ai^. Explicitly the counterterm factors 5^ t for the individual 

vertex functions read 



8, 



■WWH 

ct 



ct 



■udW,r 



5, 



■WW 

ct 



5, 



■zz 

ct 



6Z e - 




-5Z Z z + 
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sz e 
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5g 



°ct 



5. 



ZZH 




1 - 



Ml 



+ 5Z 



AZ 



+ \SZ ZZ + 8ZI - ^-5Z 



JqqZ 



Ml 



2s^ -c^5s w , 5M^ „„ 1 



r 2 



(2. 



The index r has been suppressed for those counterterms that do not depend on the chi- 
rality. The explicit expressions for the renormalization constants can, e.g., be found in 
Ref . ITS] • We merely focus on the charge renormalization constant 5Z e in the following. In 
the a(0)-scheme (i.e. the usual on-shell scheme) the electromagnetic coupling e is deduced 
from the fine-structure constant a(0), as defined in the Thomson limit. This fixes 5Z e to 



o(0) 



dZ AA (k 2 



dk 2 



% S^ z (0) 



fc 2 =0 



c w 



M| 



(2.9) 



with 



' \k 2 ) denoting the transverse part of the W gauge-boson self-energy with mo- 
mentum transfer k. In this scheme the charge renormalization constant 5Z e contains loga- 
rithms of the light-fermion masses, inducing large corrections proportional to aln(mj/s), 
which are related to the running of the electromagnetic coupling a(k 2 ) from k = to a 
high-energy scale. In order to render these quark-mass logarithms meaningful, it is nec- 
essary to adjust these masses to the asymptotic tail of the hadronic contribution to the 
vacuum polarization H AA (k 2 ) = Ti AA (k 2 )/k 2 of the photon. Using a(M|), as defined in 
Ref. [14 , as input this adjustment is implicitly incorporated, and the charge renormaliza- 
tion constant is modified to 



6Z r 



5Z t 



- -Aa(Mi) 

a(0) 2 



(2.10) 



where 



Aa(k 2 ) = U AA (0) - Re{n^ t (A; 2 )}, (2.11) 

with IIj_f t denoting the photonic vacuum polarization induced by all fermions other 
than the top quark (see also Ref. ^3]). In contrast to the a(0)-scheme the counterterm 



5Z, 



and thus the whole relative 0(a) correction in the a(M|)-scheme, does not 
involve logarithms of light quark masses, since all corrections of the form a n In™ (mj/s) are 
absorbed in the lowest-order cross section parametrized by a(Mj) = a(0)/[l — Aa(Mj)]. 
In the G M -scheme, the transition from a(0) to G M is ruled by the quantity Ar [23| IT5]. 
which is deduced from muon decay, 
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(2.12) 



Therefore, the charge renormalization constant reads 



SZ f 



5Z f 



Ar. 

«(o) 2 



(2.13) 



8 



Since Aa(Mf ) is explicitly contained in Ar, the large fermion-mass logarithms are also re- 
summed in the G^-scheme. Moreover, the lowest-order cross section in G^-parametrization 
absorbs large universal corrections to the SU(2) gauge coupling e/% induced by the p- 
parameter. 

Finally, we consider the universal corrections related to the p-parameter, or more gen- 
erally, the leading corrections induced by heavy top quarks in the loops. To this end, we 
have extracted all terms in the corrections that are enhanced by a factor m^/M^. These 
contributions are conveniently expressed in terms of 



Ap 4 



3a 



(2.14) 



which is the leading contribution to the p parameter. For the various channels we obtain 
the following correction factors to the cross sections in the G^-scheme, 



rtop 

V- 



WH 



Z0 WWH 



&qq-^ZH 



with 



&WWH 



G. 



+ A^ 



$ZZH 



2QqCw(gqqZ+ 9 qqZ ) 



qqZ) 



5 A 

•^Ap t , 
o 



(2.15) 



(2.16) 



in agreement with the results of Ref. j23]. As mentioned before, for WH production the 
only effect of Ap t is related to the WWH vertex correction in the G^-scheme, while such 
corrections to the qq'W coupling are entirely absorbed into the renormalized coupling 
e/s w . 



2.3 Real-photon emission 

Real-photonic corrections are induced by the diagrams shown in Figs. 01 and Helicity 
amplitudes for the processes qq' — > VH+j (V = W, Z) have been generated and evaluated 
using the program packages MadGraph j2H] and HELAS The result has been verified 
by an independent calculation based on standard trace techniques. The contribution a 7 
of the radiative process to the parton cross section is given by 

11 f dr 7 £ \M,\ 2 , (2.17) 

spins 



12 2s 



where the phase-space integral is defined by 

dr 7 = / To ^ / (/Z V I iS^nT WKPq+Pq'-PH-Pv - k). (2.18) 
J (2n)- i 2p Hfl J (27iy2p Vfl J (27r) J 2/c 

2.4 Treatment of soft and collinear singularities 

The phase-space integral (|2.17j) diverges in the soft (k —> 0) and collinear (p g k,Pqik — > 0) 
regions logarithmically if the photon and fermion masses are set to zero. For the treatment 
of the soft and collinear singularities we have applied the phase-space slicing method. For 
associated ZH production we have additionally applied the dipole subtraction method. 
In the following we briefly sketch these two approaches. 
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Figure 9: Bremsstrahlung diagrams for ud — > WH + 7. 




Figure 10: Bremsstrahlung diagrams for qq — ► + 7. 



2.4.1 Phase-space slicing 

Firstly, we made use of phase-space slicing, excluding the soft-photon and collinear regions 
in the integral (|2.17|) . 

In the soft-photon region m 7 < ko < AE <C the bremsstrahlung cross section 
factorizes into the lowest-order cross section and a universal eikonal factor that depends 
on the photon momentum k (see e.g. Ref. ^H])- Integration over k in the partonic CM 
frame yields a simple correction factor 5 so f t to the partonic Born cross section dcr . For 
ZH production this factor is 



8. 



soft 



ct 9 f , /2AE\ , fm 2 \ /2AE\ 1, 9 /m 2 \ fm 2 a \ 



7T - 

y 



(2.19) 



For WH production the soft factor is 

/2AE\ 



soft 



-—{Q 

2-7T I 



) 



! 

21n hd ln (^' + 21n i 



V m 7 

/2AE\ 1 , ,,/ "O 

+ o ln ^ 

V m 7 / 2 V s 



+ ln 9 



7T 



V s / 3 
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+ 



, 2AE\ ( 
2 In 1 In 



M? 



w 



, /2A£\ 1 , „ 

+ 2 In +-ln 2 , 

U 7 2 Vp^+p w 2 



Ill 



+ 2Q q 



2 In 



Pw| \(Pw + \Pw\) 2 



/2AE\ 



In 



V J 



+ Li 2 1 + 



t-M& 



(Pw ~ \Pw\) 



-2Q q 



+ Li 2 1 



2 In 



(Pw + IP 



\Y\ 



+ Li 2 ( 1 + 



In 



/2AE\ 



+ Li 2 1 + 



Pw 



IpwI) 



(2.20) 



where Q q — Q g / = +1. The difference between ZH and WH production is due to the soft 
photons emitted by the W boson. 

The factor <5 so f t can be added directly to the virtual correction factor 2 Re{<5y irt } defined 
in (|2.7|) . We have checked that the photon mass m 7 cancels in the sum 2 Re{<5^ irt } + <5 so f t . 

The remaining phase-space integration in (j2.17|) with ko > AE still contains the 
collinear singularities in the regions in which (p q k) or (pq>k) is small. Defining 0/ 7 = 
Z(p/,k) as the angle of the photon emission off / = q, q', the collinear regions are ex- 
cluded by the angular cuts 0/ 7 < AO -C 1 in the integral ()2.17|) . 

In the collinear cones the photon emission angles 0/ 7 can be integrated out. The 
resulting contribution to the bremsstrahlung cross section has the form of a convolution 
of the lowest-order cross section, 



Ccoll — C*"coll,« + CTcoll. 



tf, 



^coU,/(P/) 



2vr Jo 



dz 



111 



2z 



1 + z 2 



p ff( z )°o( z Pf)i f = q,q 

(2-21) 



with the splitting function 



(2.22) 



Note that the quark momentum p/ is reduced by the factor z so that the partonic CM 
frame for the hard scattering receives a boost. 



2.4.2 Subtraction method 

Alternatively, for ZH production, we applied the subtraction method presented in 
Ref. [2"7| . where the so-called "dipole formalism", originally introduced by Catani and 
Seymour [2E| within massless QCD, was applied to photon radiation and generalized to 
massive fermions. The general idea of a subtraction method is to subtract and to add a 
simple auxiliary function from the singular integrand. This auxiliary function has to be 
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chosen such that it cancels all singularities of the original integrand so that the phase- 
space integration of the difference can be performed numerically. Moreover, the auxiliary 
function has to be simple enough so that it can be integrated over the singular regions 
analytically, when the subtracted contribution is added again. 

The dipole subtraction function consists of contributions labelled by all ordered pairs 
of charged external particles, one of which is called emitter, the other one spectator. For 
qq — > ZH we, thus, have 2 different emitter/spectator cases //': qq, qq. The subtraction 
function that is subtracted from Sspins I M--* | 2 is given by 



\M 



sub | 



Q 2 q e 2 \g q q(p qj p q ,k) \Ml(x qq p q ,Pq,kz. 



VI) 



spins 



+ 9qq(Pq, Pq, J2 \ M o(Pq, X qqPq, k Z. 



qq) 



(2.23) 



spins 



with the functions 

9qq(Pq,Pq, 
9qq(Pq,Pq,k) 

and the auxiliary variable 



1 



(p q k)x qq 
1 



1 — x 



X, 



qq 



1 — X, 



X 



qq 



(Pqk)x qq 

PqPq ~ Pqk ~ Pqk 
PqPq 



X 



qq 



qq 



(2.24) 



(2.25) 



For the evaluation of |.M su b| 2 in (|2.23|) the Z-boson momenta kzjf still have to be speci- 
fied. They are given by 



K Z,qq 



A(p q , Pq Y v k" z , 
with the Lorentz transformation matrix 



K z, qq 



A{pq,p q T v k 



(2.26) 



A(pi,P2) 



9' 



(P + PY(P + P)„ 2P»P V 



P 2 + PP 



p2 



P» = k£ + k% 1 , P» = x qqP Z+ P %. 

(2.27) 

The modified Z-boson momenta kzjf still obey the on-shell condition k\ f*, = M|, and 
the same is true for the corresponding Higgs-boson momenta that result from momentum 
conservation. It is straightforward to check that all collinear and soft singularities cancel 
in Lupins I-M7I 2 — |.M su b| 2 so that this difference can be integrated numerically over the 
entire phase space (j2.18|) . 

The contribution of |A^ su b| 2 , which has been subtracted by hand, has to be added 
again. This is done after the singular degrees of freedom in the phase space (j2.18|) are 
integrated out analytically, keeping an infinitesimal photon mass m 7 and small fermion 
masses mj as regulators j2jj. The resulting contribution is split into two parts: one that 
factorizes from the lowest-order cross section &o and another part that has the form of a 
convolution integral over <5"o with reduced CM energy. The first part is given by 



Osub,l 



Q 2 — 



2£(s,m 2 ) + 3 



2vr 2 



00 



(2.28) 
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with the auxiliary function 



(mi 



V r 



(mi 



V r 



(mi 



C(r,rrii) = hi In — Y - + ln — M --In 2 -£ + - In -i . 



V r 



V r 



(mi 



V r 



(2.29) 



The IR and fermion-mass singularities contained in d<5" su b,i exactly cancel those of the 
virtual corrections. The second integrated subtraction contribution is given by 



^sub,2(Pg,Pg) = Q\l^ j dx {[Gqq{s , X ))+ &o{ X Pq, Pq) + \Qq q (§, x)) + a {p q , Xp q ) } 



+ 



Q 2 q7T I dx(l-x)fcr (xp q ,pg) +(ro{p q ,xpg 

Z7T JO I T q^- T q 



Ta—>-Tn 



where the usual 



prescription, 



da; 



/(*) 



dxf(x) \g(x) - g(l)] , 



is applied to the integration kernels 

Gqq(r,x) = Gqq(r,x) = P ff (x) 



in -L 



S=(p 9 +Pg) 2 



(2.30) 



(2.31) 



(2.32) 



In ()2.3())1 we have indicated explicitly how the Mandelstam variable r has to be chosen in 
terms of the momenta in the evaluation of the part containing [Gff>(r, x)] + . Note, however, 
that in ()2.30|) the variable s that is implicitly used in the calculation of <5"o(. . .) is reduced 

tO 2xp q Pq = XS. 

In summary, within the subtraction approach the real correction reads 



11 / dr. 

12 2S 



\M 



sub | 



spins 



+ C"sub,l + C"sub,2- 



(2.33) 



It should be realized that in <7 su b,i and <3" su b,2 the full photonic phase space is integrated 
over. This does, however, not restrict the subtraction approach to observables that are 
fully inclusive with respect to emitted photons, but rather to observables that are inclusive 
with respect to photons that are soft or collinear to any charged external fermion (see 
discussions in Sect. 6.2 of Ref. [2Zj and Sect. 7 of Ref. 



3 The hadron cross section 

The proton-(anti-)proton cross section a is obtained from the parton cross sections a^ qiq2 ^ 
by convolution with the corresponding parton distribution functions qi^{x), 

da(s) = V I dxi / dx2qi(x 1 )q 2 (x2)dcr {qiq2 \pq 1 ,pq 2 ), (3.1) 
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where X\^ are the respective momentum fractions carried by the partons q\^- In the sum 
Y, qiq2 the quark pairs qiq 2 run over all possible combinations qq' and q'q where q = u, c and 
q' = d,s for WH production and q = q' = u, d, s, c, b for ZH production. The squared CM 
energy s of the pp (pp) system is related to the squared parton CM energy s by s = x\x 2 s. 

The (9(a)-corrected parton cross section a^ 1 ^ contains mass singularities of the form 
aln(ra g ), which are due to collinear photon radiation off the initial-state quarks. In com- 
plete analogy to the MS factorization scheme for next-to-leading order QCD corrections, 
we absorb these collinear singularities into the quark distributions. This is achieved by 
replacing q(x) in (|3.1|) according to 



q(x) — > q(x, M 2 



1 dz fx , ,« 
— q(-,M 2 

x Z \Z 



a ^ 2 f / M 2 
— Q\ In — 
Z7r L V m z q 



P 



ff\ 



P„(z) 21n(l-z) + l 



(3.2) 



where M is the factorization scale (see Ref. llj). This replacement defines the same finite 
parts in the O(o) correction as the usual MS factorization in Z)-dimensional regularization 
for exactly massless partons, where the ln(m ? ) terms appear as 1/(D — 4) poles. In (J3.2)) 
we have regularized the soft-photon pole by using the [...]+ prescription. This procedure 
is fully equivalent to the application of a soft-photon cutoff AE (see Ref. [TU]) where 



q(x) — > q(x, M 2 



- - Q 2 Al - ln(2A J E/v / I) - ln 2 (2AE/v / I) 

7r y 

ln(2AB/-v/J) + |) IngJ)} 



l-2AE/Vldz fx a ^ 2 „ . . f, (M 2 



(1 



(3.3) 



The absorption of the collinear singularities of 0(a) into quark distributions, as a 
matter of fact, requires also the inclusion of the corresponding 0(a) corrections into the 
DGLAP evolution of these distributions and into their fit to experimental data. At present, 
this full incorporation of 0(a) effects in the determination of the quark distributions has 
not yet been performed. However, an approximate inclusion of the 0(a) corrections 
to the DGLAP evolution shows that the impact of these corrections on the quark 
distributions in the MS factorization scheme is well below 1%, at least in the x range 
that is relevant for associated VH production at the Tevatron and the LHC. Therefore, 
the neglect of these corrections to the parton distributions is justified for the following 
numerical study. 
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4 Numerical results 



4.1 Input parameters 

For the numerical evaluation we used the following set of input parameters [30j, 



Gfj, 


= 1.16639 x lO^GeV -2 , 


a(0) 


= 1/137.03599976, 


a{M z ) 


= 1/128.930, 


M w 


= 80.423 GeV, 


M z 


= 91.1876 GeV, 






m e 


= 0.510998902 MeV, 


«V 


= 105.658357 MeV, 


m T 


= 1.77699 GeV, 


m a 


= 66 MeV, 


m c 


= 1.2 GeV, 


m t 


= 174.3 GeV, 


m d 


= 66 MeV, 


m s 


= 150 MeV, 


m h 


= 4.3 GeV 


|Kd| 


= 0.975, 


\v w \ 


= 0.222, 






|Kd| 


= 0.222, 


\v CB \ 


= 0.975. 







(4.1) 

The masses of the light quarks are adjusted such as to reproduce the hadronic contribution 
to the photonic vacuum polarization of Ref . Jl] . They are relevant only for the evaluation 
of the charge renormalization constant 5Z e in the a(0)-scheme. For the calculation of 
the pp and pp cross sections we have adopted the CTEQ6L1 and CTEQ6M [HI] parton 
distribution functions at LO and 0(a s ), corresponding to A^ = 165 MeV and Ag = 
226 MeV at the one- and two-loop level of the strong coupling a s (fi), respectively. The top 
quark is decoupled from the running of a s (n). If not stated otherwise the factorization 
scale M is set to the invariant mass of the Higgs-vector-boson pair, M = y/syn- For the 
treatment of the soft and collinear singularities we have applied the phase-space slicing 
method as described in Sect. 12.41 We have verified that the results are independent of 
the slicing parameters 2AE/\^I and A9 when these parameters are varied within the 
range 10~ 2 — 10~ 4 . In the case of associated ZH production we have, in addition, applied 
the dipole subtraction method. The results agree with those obtained using phase-space 
slicing. We observe that the integration error of the subtraction method is smaller than 
that of the slicing method by at least a factor of two. 

4.2 Electroweak corrections 

In this subsection we present the impact of the electroweak 0(a) corrections on the cross 
section predictions for the processes pp/pp — > W + H + X and pp/pp — > ZH + X at the 
Tevatron and the LHC Figures ITD and fT2l show the relative size of the O(a) corrections 
as a function of the Higgs-boson mass for pp — > W + H + X and pp — > ZH + X at the 
Tevatron. Results are presented for the three different input-parameter schemes. The 
corrections in the G^- and «(M|)-schemes are significant and reduce the cross section by 
5-9% and by 10-15%, respectively. The corrections in the a(0)-scheme differ from those in 
the G M -scheme by 2Ar w 6% and from those in the a(M|)-scheme by 2Aa(M|) w 12%. 
The fact that the relative corrections in the a(0)-scheme are rather small results from 
accidental cancellations between the running of the electromagnetic coupling, which leads 
to a contribution of about 2Aa(M|) w +12%, and other (negative) corrections of non- 
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universal origin. Thus, corrections beyond 0{a) in the a(0)-scheme cannot be expected 
to be suppressed as well. In all schemes, the size of the corrections does not depend 
strongly on the Higgs-boson mass. The unphysical singularities at the thresholds M# = 
2Mw and 2Mz can be removed by taking into account the finite widths of the unstable 
particles, see e.g. Refs. Representative results for the leading-order cross section and 
the electroweak 0(a) corrections are collected in Tables [T] and EJ 

Figures IT31 and ITU and Tables 01 and 0] show the corresponding results for pp —>■ W + H + 
X and pp — > ZH + X at the LHC. The corrections are similar in size to those at the 
Tevatron and reduce the cross section by 5-10% in the G^-scheme and by 12-17% in the 
a(M|)-scheme. We note that the electroweak corrections to pp — > W~ H + X at the LHC 
differ from those to pp — > W + H + X by less than about 2%. 

In order to unravel the origin of the electroweak corrections we display the contri- 
butions of individual gauge-invariant building blocks. Figure separates the fermionic 
corrections (comprising all diagrams with closed fermion loops) from the remaining bosonic 
contributions to pp —>■ W + H + X at the Tevatron in the G M -scheme. We observe that 
the bosonic corrections are dominant and that bosonic and fermionic contributions partly 
compensate each other. A similar result is found for the pp — ► ZH+X cross section, where 
we display the gauge-invariant contributions from (photonic) QED corrections, fermionic 
corrections, and weak bosonic corrections in Fig. Note that large logarithmic correc- 
tions from initial-state photon radiation have been absorbed into the quark distribution 
functions. The remainder of the QED corrections turns out to be strongly suppressed with 
respect to the fermionic and weak bosonic corrections. A similar pattern is observed for 
the pp — > W + H + X and pp — > ZH + X cross sections at the LHC, see Figs. ITTlandfTKl 

At first sight, the large size of the non-universal corrections, i.e. corrections that are 
not due to the running of a(k 2 ), photon radiation, or other universal effects, might be 
surprising. However, a similar pattern has already been observed in the electroweak cor- 
rections to the processes e + e~ — ► ZH ^3] and e + e~ — > Z*H — > vvH p!3lll6| . Also there 
large non-universal fermionic and bosonic corrections of opposite sign occur. It was also 
observed that the corrections cannot be approximated by simple formulae resulting from 
appropriate asymptotic limits. 2 For instance, taking the large top-mass limit (m t — ■> oo) 
in the fermionic corrections to WH production in the G^-scheme (see Sect. I2.2~2*)) . the 



~ —1.6%, which even differs 



leading term in the relative correction is given by S q ^_+ WH 

in sign from the full result (see Fig. I15j) . The reason for this failure is that the relevant 
scale in the WWH vertex, from which the leading m 2 term in the limit m t —> oo results, 
is set by the variable s which is not much smaller than but rather of the same order as 
2 . For the ZH channel, we get 8°?_+ ZH ~ ^TI-^zh ~ —1%, again reflecting the 



mt 



G. 



Gri 



failure of the heavy-top limit as a suitable approximation. Concerning the weak bosonic 
corrections, large negative contributions are expected in the high-energy limit owing to the 
occurrence of Sudakov logarithms of the form — a/7rlog 2 (s/M 2 v ). However, the relevant 

2 In Ref. [31] the part of the fermion- loop correction that is enhanced by an explicit factor am^/M^f 
was calculated. Moreover, in the second paper of Ref. also diagrams with internal Higgs bosons were 
taken into account. Using a = 1/128, which roughly corresponds to the a(Af|)-scheme, these authors find 
about —1% to —2% for the sum of these corrections, which were assumed to be the leading ones. This 
has to be compared with our result of about —12% for the full O(a) corrections in the a(M|)-scheme. 
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partonic CM energies are not yet large enough for the Sudakov logarithms to provide 
a good approximation for the full corrections. 

4.3 The cross section at NLO 

In this subsection we present the cross section prediction for associated WH and ZH 
production at the Tevatron and at the LHC, including the NLO order electroweak and 
QCD corrections, and we quantify the residual theoretical uncertainty due to scale vari- 
ation and the parton distribution functions. The total cross sections for the processes 
pp/pp — > W ± H + X (sum of W + H and W~H) and pp/pp — > ZH + X at the Tevatron 
and the LHC are displayed in Figs. [TI5H221 Representative results are listed in Tables [SHE 
For the central renormalization and factorization scale fj, = = ^Jsvh the NLO QCD 
corrections increase the LO cross section by typically 20-25%. As discussed in detail in 
Sect. 14.21 the NLO electroweak corrections are sizeable and decrease the cross section by 
5-10% in the G^-scheme. The size of the 0(a B ) and 0(a) corrections does not depend 
strongly on the Higgs-boson mass. 

The NLO prediction is very stable under variation of the QCD renormalization and 
factorization scales. We have varied both scales independently in the range //o/5 < \i < 
5/io- For both the Tevatron and the LHC, the cross section increases monotonically with 
decreasing renormalization scale. At the Tevatron, the maximal (minimal) cross section is 
obtained choosing both the renormalization and factorization scales small (large). At the 
LHC, in contrast, the maximal (minimal) cross section corresponds to choosing a large 
(small) factorization scale. From the numbers listed in Tables EHH1 one can conclude that 
the theoretical uncertainty introduced by varying the QCD scales in the range /io/5 < [i < 
5/io is less than approximately 10%. We have verified that the QED factorization-scale 
dependence of the C(a)-corrected cross section is below 1% and thus negligible compared 
to the other theoretical uncertainties. The QED scale dependence should be reduced 
further when using QED-improved parton densities. 

We have also studied the uncertainty in the cross-section prediction due to the error 
in the parametrization of the parton densities. To this end we have compared the NLO 
cross section evaluated using the default CTEQ6 [3T] parametrization with the cross sec- 
tion evaluated using the MRST2001 35 j parametrization. The results are collected in 
Tables IHHT21 Both the CTEQ and MRST parametrizations include parton-distribution- 
error packages which provide a quantitative estimate of the corresponding uncertainties 
in the cross sections. 3 Using the parton-distribution-error packages and comparing the 
CTEQ and MRST2001 parametrizations, we find that the uncertainty in predicting the 
processes pp/pp — > W ± H + X and pp/pp — > ZH + X at the Tevatron and the LHC due 
to the parametrization of the parton densities is less than approximately 5%. 

3 In addition, the MRST (21) parametrization allows to study the uncertainty of the NLO cross section 
due to the variation of a s . For associated WH and ZH hadroproduction, the sensitivity of the theoretical 
prediction to the variation of a s (a s (M|) = 0.119 ± 0.02) turns out to be below 2%. 
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5 Conclusions 



We have calculated the electroweak 0(a) corrections to Higgs-boson production in associ- 
ation with W or Z bosons at hadron colliders. These corrections decrease the theoretical 
prediction by up to 5-10%, depending in detail on the Higgs-boson mass and the input- 
parameter scheme. We have updated the cross section prediction for associated WH 
and ZH production at the Tevatron and at the LHC, including the next-to-leading order 
electroweak and QCD corrections. Finally, the remaining theoretical uncertainty has been 
studied by varying the renormalization and factorization scales and by taking into account 
the uncertainties in the parton distribution functions. We find that the scale dependence 
is reduced to about 10% at next-to-leading order, while the uncertainty due to the parton 
densities is less than about 5%. 
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Figure 11: Relative electroweak correction S as a function of the Higgs-boson mass for the 
total cross section pp — > W + H + X (-^/s = 1.96 TeV). Results are presented for the a(0)-, 
a(M|)-, and G M - schemes. 
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Figure 12: Relative electroweak correction S as a function of the Higgs-boson mass for the 
total cross section pp — > ZH + X (y/s = 1.96 TeV). 
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Figure 13: Relative electroweak correction 5 as a function of the Higgs-boson mass for 
the total cross section pp — > W + H + X {yfs = 14 TeV). 
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Figure 14: Relative electroweak correction 5 as a function of the Higgs-boson mass for the 
total cross section pp — >■ ZH + X {-sfs = 14 TeV). 
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Figure 15: Different contributions to the relative electroweak correction 5 in the G M - 
scheme as a function of the Higgs-boson mass for the total cross section pp — > W + H + X 
{yfs = 1.96 TeV). 
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Figure 16: Different contributions to the relative electroweak correction 5 in the G M - 
scheme as a function of the Higgs-boson mass for the total cross section pp — > ZH + X 
(y/s = 1.96 TeV). 
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Figure 17: Different contributions to the relative electroweak correction 8 in the G M - 
scheme as a function of the Higgs-boson mass for the total cross section pp — > W + H + X 
(y/s = 14 TeV). 
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Figure 18: Different contributions to the relative electroweak correction 5 in the G M - 
scheme as a function of the Higgs-boson mass for the total cross section pp — > ZH + X 
(y/s = 14 TeV). 
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Figure 19: Total cross section for pp — ► W^H + X (sum of W + H and W~H) at the 
Tevatron (y/s = 1.96 TeV) in LO, NLO QCD, and including NLO QCD and electroweak 
corrections in the G^-scheme. The renormalization and factorization scales have been set 
to the invariant mass of the Higgs-vector-boson pair, fi = y/syu- CTEQ6L1 and CTEQ6M 
pTT] parton distribution functions have been adopted at LO and 0(a s ), respectively. 
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Figure 20: Total cross section for pp — > ZH + X at the Tevatron (y/s = 1.96 TeV) in LO, 
NLO QCD, and including NLO QCD and electroweak corrections in the G^-scheme. 
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Figure 21: Total cross section for pp -> W ± H + X at the LHC (y/s = 14 TeV) in LO, 
NLO QCD, and including NLO QCD and electroweak corrections in the G^-scheme. 



pp ->• ZH + X 

= 14 TeV 




M„ / GeV 

rl 

Figure 22: Total cross section for pp -> ZH + X at the LHC (y/s = 14 TeV) in LO, NLO 
QCD, and including NLO QCD and electroweak corrections in the G^-scheme. 
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140.00 


0.02945(1) 


-0.89(1) 


0.03327(1) 


-13.53(1) 


0.03149(1) 


-7.05(1) 


170.00 


0.01367(1) 


-2.55(1) 


0.01544(1) 


-15.30(1) 


0.01461(1) 


-8.90(1) 


190.00 


0.008517(1) 


-1.97(1) 


0.009624(1) 


-14.66(1) 


0.009106(1) 


-8.38(1) 



Table 1: Total lowest-order hadronic cross section a (pp — * W + H+X) and corresponding 
relative electroweak correction 5 (y/s = 1.96 TeV). Results are presented for the ct(0)-, 
a(M|)-, and G^-schemes. The integration error is given in brackets. 



M H / GeV 


<7o|a(o)/pb 


5\a(0)/% 


o"oU(M|)/pb 




ff o|c M /pb 


S\gJ% 


80.00 


0.2199(1) 


0.99(1) 


0.2484(1) 


-11.52(1) 


0.2350(1) 


-4.73(1) 


100.00 


0.1142(1) 


0.95(1) 


0.1290(1) 


-11.57(1) 


0.1221(1) 


-4.91(1) 


120.00 


0.06358(1) 


0.97(1) 


0.07182(1) 


-11.55(1) 


0.06796(1) 


-5.01(1) 


140.00 


0.03727(1) 


0.91(1) 


0.04211(1) 


-11.61(1) 


0.03984(1) 


-5.18(1) 


170.00 


0.01799(1) 


1.12(1) 


0.02032(1) 


-11.38(1) 


0.01922(1) 


-5.10(1) 


190.00 


0.01148(1) 


1.26(1) 


0.01297(1) 


-11.24(1) 


0.01227(1) 


-5.04(1) 



Table 2: Total lowest-order hadronic cross section ao(pp — > Z H + X) and corresponding 
relative electroweak correction 5 (y/s = 1.96 TeV). 



27 



Mh 1 Ge V 


0"o|a(o)/pb 


5\a(0)/% 


°"0 U(Af|)/P b 


<*la(M§)/% 


ff o|G M /pb 


S\gJ% 


80.00 


2.660(1) 


0.31(1) 


3.005(1) 


-12.22(2) 


2.844(1) 


-5.43(1) 


100.00 


1.410(1) 


-0.11(1) 


1.594(1) 


-12.67(2) 


1.508(1) 


-5.99(1) 


120.00 


0.8114(2) 


-0.65(1) 


0.9166(2) 


-13.24(2) 


0.8673(2) 


-6.67(1) 


140.00 


0.4967(1) 


-1.49(1) 


0.5610(1) 


-14.16(2) 


0.5309(1) 


-7.68(1) 


170.00 


0.2605(1) 


-3.33(1) 


0.2942(1) 


-16.12(2) 


0.2784(1) 


-9.72(2) 


190.00 


0.1776(1) 


-2.92(1) 


0.2007(1) 


-15.67(2) 


0.1899(1) 


-9.36(1) 



Table 3: Total lowest-order hadronic cross section <7o{pp —* W + H+X) and corresponding 
relative electroweak correction 5 (y/s = 14 TeV). 



M H / Ge V 


°0 a(0)/pb 


5\a(0)/% 


°"oU(M|)/pb 


^U(M|)/% 


ff o|G M /pb 


S\gJ% 


80.00 


2.299(1) 


0.95(1) 


2.595(1) 


-11.56(1) 


2.457(1) 


-4.77(1) 


100.00 


1.232(1) 


0.83(1) 


1.392(1) 


-11.68(1) 


1.317(1) 


-5.03(1) 


120.00 


0.7134(1) 


0.76(1) 


0.8058(1) 


-11.77(1) 


0.7630(1) 


-5.22(1) 


140.00 


0.4381(1) 


0.54(1) 


0.4950(1) 


-12.01(1) 


0.4684(1) 


-5.56(1) 


170.00 


0.2297(1) 


0.37(1) 


0.2595(1) 


-12.18(1) 


0.2456(1) 


-5.88(1) 


190.00 


0.1563(1) 


0.32(1) 


0.1765(1) 


-12.23(1) 


0.1670(1) 


-6.01(1) 



Table 4: Total lowest-order hadronic cross section cr (pp — > Z H + X) and corresponding 
relative electroweak correction 5 (y/s = 14 TeV). 
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Mu/GeV 

J, V-L J-l I VJ \j V 


(Tn n / nil 


u NLO / t> u 


aQCD+EW/ b 


QCD+EW / i 
^NLO /P b 


= 5/i 
/j F = 5/i 


Vr = A*o/5 
= /V 5 


oU.UU 


n /i 1 1 7/1 ^ 
U.411 / (1 ) 




11.04114^/ J 


0.5033(1) 


0.5838(1) 


100.00 


0.2056(1) 


0.2801(1) 


0.2685(1) 


0.2482(1) 


0.2911(1) 


120.00 


0.1106(1) 


0.1504(1) 


0.1436(1) 


0.1318(1) 


0.1562(1) 


140.00 


0.06297(1) 


0.08536(1) 


0.08092(1) 


0.07377(1) 


0.08833(1) 


170.00 


0.02921(1) 


0.03940(1) 


0.03679(1) 


0.03318(1) 


0.04037(1) 


190.00 


0.01821(1) 


0.02446(1) 


0.02294(1) 


0.02056(1) 


0.02525(1) 



Table 5: Total cross section for pp —>■ W ± H + X (sum of W + H and W~H) at the Tevatron 
(y/s = 1.96 TeV) in LO, NLO QCD, and including NLO QCD and electroweak corrections 
in the G^-scheme. The renormalization scale (jir) and the factorization scale (jif) have 
been set to the invariant mass of the Higgs-vector-boson pair, /i = /i = y/syu- CTEQ6L1 
and CTEQ6M [HI] parton distribution functions have been adopted at LO and 0(a s ), 
respectively. The last two columns show the minimal and maximal cross section prediction 
obtained from varying the QCD renormalization and factorization scales independently in 
the range /io/5 < \i < 5/io- 



M H /GeV 


ff o|G M /pb 


QCD / , 
ff NLO /Pb 


QCD+EW / i 
^NLO /Pb 


QCD+EW / i 
^NLO /Pb 


Vr — 5/io 
/i F = 5/i 


Hr — /^o/S 
= /io/5 


80.00 


0.2350(1) 


0.3181(1) 


0.3070(1) 


0.2858(1) 


0.3317(1) 


100.00 


0.1221(1) 


0.1649(1) 


0.1589(1) 


0.1470(1) 


0.1722(1) 


120.00 


0.06796(1) 


0.09160(1) 


0.08820(1) 


0.08111(2) 


0.09575(1) 


140.00 


0.03984(1) 


0.05354(1) 


0.05148(1) 


0.04706(2) 


0.05604(1) 


170.00 


0.01922(1) 


0.02570(1) 


0.02472(1) 


0.02242(1) 


0.02701(1) 


190.00 


0.01227(1) 


0.01635(1) 


0.01573(1) 


0.01418(1) 


0.01722(1) 



Table 6: Total cross section for pp — > ZH + X at the Tevatron (y/s = 1.96 TeV) in LO, 
NLO QCD, and including NLO QCD and electroweak corrections in the G^-scheme. 
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Mrr /npV 

lvljjl Vjrc V 




rr QCD /nh 
a NLO/P D 


QCD+EW , h 
a NLO /P D 


QCD+EW / K 
°NLO /P b 


HR = 5yU 

Mf = Mo/5 


Mi? = Mo/5 
Mf = 5/i 


80.00 


4.679(2) 


5.676(2) 


5.423(2) 


4.875(2) 


5.749(5) 


100.00 


2.462(1) 


3.005(1) 


2.859(1) 


2.606(2) 


3.033(2) 


120.00 


1.405(1) 


1.726(1) 


1.633(1) 


1.505(1) 


1.731(1) 


140.00 


0.8537(2) 


1.054(1) 


0.9892(3) 


0.9204(3) 


1.050(1) 


170.00 


0.4434(1) 


0.5504(1) 


0.5078(1) 


0.4782(1) 


0.5388(3) 


190.00 


0.3003(1) 


0.3745(1) 


0.3466(1) 


0.3285(1) 


0.3675(2) 



Table 7: Total cross section for pp -> W^if + X at the LHC (y/s = 14 TeV) in LO, NLO 
QCD, and including NLO QCD and electroweak corrections in the G^-scheme. 



M H /GeV 


^o|c M /pb 


QCD / , 
^NLO / P b 


QCD+EW / i 
^NLO /Pb 


QCD+EW / i 

°"nlo /Pb 


M-R — 5/io 
Mf = Mo/5 


Mr — Mo/5 
Mf = 5mo 


80.00 


2.457(1) 


2.974(1) 


2.857(1) 


2.578(2) 


3.018(3) 


100.00 


1.317(1) 


1.605(1) 


1.539(1) 


1.407(1) 


1.629(1) 


120.00 


0.7630(1) 


0.9346(3) 


0.8947(3) 


0.8271(2) 


0.9462(6) 


140.00 


0.4684(2) 


0.5768(2) 


0.5508(2) 


0.5138(1) 


0.5830(3) 


170.00 


0.2456(1) 


0.3045(1) 


0.2900(1) 


0.2736(1) 


0.3068(2) 


190.00 


0.1670(1) 


0.2078(1) 


0.1978(1) 


0.1879(1) 


0.2094(1) 



Table 8: Total cross section for pp -> ZH + X at the LHC (^s = 14 TeV) in LO, NLO 
QCD, and including NLO QCD and electroweak corrections in the G^-scheme. 
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M H /GeV 


CTEQ6M [H] 


MRST2001 [SS] 


80.00 


0.5404(2) ± 0.021 


0.5448(2) ± 0.0097 


100.00 


0.2685(1) ± 0.011 


0.2698(1) ± 0.0052 


120.00 


0.1436(1) ± 0.0060 


0.1437(1) ± 0.0030 


140.00 


0.08092(1) ± 0.0035 


0.08065(1) ± 0.0018 


170.00 


0.03679(1) ± 0.0017 


0.03644(1) ± 0.00091 


190.00 


0.02294(1) ± 0.0011 


0.02262(1) ± 0.00060 



Table 9: Parton distribution function (PDF) uncertainties: Total cross section for 
pp — ► W ± H + X at the Tevatron (*J~s = 1.96 TeV) including NLO QCD and electroweak 
corrections in the G^-scheme for different sets of parton distribution functions. The results 
include an estimate of the uncertainty due to the parametrization of the parton densities 
as obtained with the CTEQ6 [3T] and MRST2001 jSHj eigenvector sets (columns 2 and 
3, respectively). The renormalization and the factorization scales have been set to the 
invariant mass of the Higgs-vector-boson pair, // = /i = y/svH- 



M H /GeV 


CTEQ6M |H] 


MRST2001 |35] 


80.00 


0.3070(1) ± 0.012 


0.3090(1) ± 0.0039 


100.00 


0.1589(1) ± 0.0064 


0.1596(1) ± 0.0020 


120.00 


0.08820(1) ± 0.0036 


0.08840(1) ± 0.0011 


140.00 


0.05148(1) ± 0.0021 


0.05151(1) ± 0.00066 


170.00 


0.02472(1) ± 0.0010 


0.02469(1) ± 0.00033 


190.00 


0.01573(1) ± 0.00068 


0.01568(1) ± 0.00021 



Table 10: PDF uncertainties: Total cross section for pp — > ZH + X at the Tevatron 
(y/s = 1.96 TeV) including NLO QCD and electroweak corrections in the G M -scheme for 
different sets of parton distribution functions. 
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Mu 1 Ge V 


CTEQ6M [HI] 


MRST2001 


80.00 


5.423(2) ± 0.18 


5.509(2) ± 0.071 


100.00 


2.859(1) ± 0.096 


2.910(1) ± 0.035 


120.00 


1.633(1) ± 0.055 


1.664(1) ± 0.021 


140.00 


0.9892(3) ± 0.034 


1.010(1) ± 0.012 


170.00 


0.5078(1) ± 0.018 


0.5193(1) ± 0.0063 


190.00 


0.3466(1) ± 0.012 


0.3547(2) ± 0.0043 



Table 11: PDF uncertainties: Total cross section for pp — > W ± H + X at the LHC (y/s = 
14 TeV) including NLO QCD and electroweak corrections in the G^-scheme for different 
sets of parton distribution functions. 



M H /GeV 


CTEQ6M |SI] 


MRST2001 [SS! 


80.00 


2.857(1) ± 0.095 


2.936(1) ± 0.036 


100.00 


1.539(1) ± 0.051 


1.583(1) ± 0.019 


120.00 


0.8947(3) ± 0.030 


0.9217(3) ± 0.011 


140.00 


0.5508(2) ± 0.019 


0.5681(2) ± 0.0067 


170.00 


0.2900(1) ± 0.010 


0.2994(1) ± 0.0036 


190.00 


0.1978(1) ± 0.0069 


0.2045(1) ± 0.0025 



Table 12: PDF uncertainties: Total cross section for pp — > ZH + X at the LHC (\/s = 
14 TeV) including NLO QCD and electroweak corrections in the G^-scheme for different 
sets of parton distribution functions. 
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